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(— I Abstract 

I Quantum Hall systems are a suitable theme for a case study in the general area of nan- 

D otechnology. In particular, it is a good framework in which to search for universal principles 

relevant to nanosystem modeling, and nanosystem-specific signal processing. Recently, we 

O have been able to construct a PDE model of a quantum Hall system, which consists of the 

^ Schrodinger equation supplemented with a special type nonlinear feedback loop. This result 

^S^ stems from a novel theoretical approach, which in particular brings to the fore the notion of 

1 ^^ quantum information. In this article we undertake to modify the original model by substi- 

—^ tuting the dynamics based on the Dirac operator. This leads to a model that consists of a 

<*' system of three nonlinear ly coupled first order elliptic equations in the plane. 

OO 

cn 

\o 

^. 1 Quantum entanglement of composite systems and its 

"^ 

§ associated Hamiltonian dynamics 

o 

|> In this section we give a very brief overview of the principles behind the mesoscopic loop models 
^ as developed in [H El |3l IH [5] . The theme of quantum information is extremely relevant to our 
c3 approach, but it has been brought to light in explicit terms only recently; an in-depth discussion 
will appear elsewhere, [6j. 

Since the measurement of Hall resistance in a quantum-Hall system results in a classical signal, 
it is natural to ask if any quantum information is in the process transduced to the classical medium. 
Let us develop this point of view into a formal discussion. Suppose that the electronic solid-state 
system and the ambient magnetic field are represented by two Hilbert spaces, say, H and H 
respectively. According to the principles of quantum mechanics, [7], the combined system is then 



represented by the tensor product space H H. Moreover, the states of the composite system 
assume the form 

l^combined) = Y.K^n\^m) ® |^n) G H(g)H. (1) 

We will associate with every state as above an entanglement operator (with no intention to imply 
that the composite state ([I]) is necessarily entangled in the rigorous sense of the term, [7]) 

K = ^AmnWm){iJn\- (2) 

Thus, if an electronic system is found in the state V" ^ H, then the electromagnetic system will 
collapse onto the state 

Let us define the density operator of system H in the usual way 

P ~ fj V combined'' combined I J ■ ^ ' 

One readily finds 

p = K*K. (4) 

This lays out the vocabulary for a quantum-mechanical description of a composite system. Next, 
we postulate a model relevant specifically to the quantum Hall systems. Namely, we propose that 
in the presence of a magnetic field, the dynamics of the electronic system is derived from the 
energy functional of the form 

Tr (ifp)+/? log det(p), 

where H is the "regular" Hamiltonian governing the microscopic dynamics of the electronic sub- 
system. For the purposes at hand there is no need to discuss the constant /3 beyond the simple 
statement that it is nonzero when the magnetic field is nonvanishing. Furthermore, we assume for 
simplicity that H is diagonalizable in the basis consisting of its eigenvalues 

Hi)k = Eki)k- (5) 

Now, substituting the decomposition Q into the above yields the following Hamiltonian 

Z{K) = Tr (KHK*) + (3 log det {K*K) . (6) 

This in turn leads to the following system of Hamiltonian equations, see [T] and [2]: 

-ihk = KH + (3K*-\ (7) 

(From this point on we will set h = 1.) The corresponding eigenvalue problem assumes the form 

KH + (3K*-^ = uK. (8) 



We postulate that homeostasis of the quantum Hall system is captured by the latter equation. 
Solutions of ([s]) are found to assume the form (P, [2]): 

k=y: (-^Y \u^,){i^,\. (9) 

where [/ : H — >^ H is an arbitrary unitary transformation. 

2 The mesoscopic-loop model based on the Dirac-type 
equations 

The model we are about to present deals with two types of particles, which will be called electrons 
and holes. Both of them will have the same effective mass m* = 1, and charge e = 1. It seems 
appropriate to emphasize that, of course, the notion of effective mass depends on the model. 
All analysis will be performed in a two-dimensional plane. Recall that the Hodge star * in two 
dimensions is a linear operation on differential forms determined by the relations *dx = dy, *dy = 
—dx, *1 = dx A dy, *{dx A dy) = 1. The exterior derivative d defines the co-derivative 6 = — * d*. 
Let the vector potential A be given in local coordinates as A = Aidx + A2dy. In particular, 
6A = -(^1,^ + ^2 J, while 

B{x,y) = *dA = A2,,-A,^y. (10) 

is the magnetic flux density. The correspondng Schrodinger operator 5*^ acts on a wavefunction 
ip as follows 

25a^ = -{dl + dl)^ + 2z{Aid,ij + A2dy^) + (A? + A^)^ - z(M)^ (11) 

(The factor of 2 in front of Sa stems from the fact that in the adopted units the kinetic energy is 
multiplied by a factor of fi^ /2m* = 1/2.) The Schrodinger operator can be successfully used in a 
construction of mesoscopic loop models, [S] . In this article we will consider a modification based 
on a Dirac type factorization of the Schrodinger operator. To this end we use the two-dimensional 
Dirac operator, [8J, with its parameters set at m = and c = 1. Namely, let 



Ha 

where 



D* 
D 



(12) 



D = -id, -A^ + dy- iA2, (13) 

so that 

D* = -id, - Ai - dy + iA2. (14) 



A direct calculation shows that 



Hi 



D*D 
DD* 



2Sa-B 

2Sa + B 



(15) 



where B = B(x,y) is determined by A via (10). Observe the difference of signs at the magnetic 



flux density term. This necessitates the interpretation of the two wavefunctions defining the 
statevector 

as representing particles with opposite pseudo-spins. At this point we venture to describe a 
quantum Hall system in a manner discussed in the previous section, which brings the notion of 
entanglement to the fore. Namely, we consider the following system of equations 



H, 



i^2 



i^dA{x,y) 



E 



\K 



^1 ' 

^2 / 

^2 



(16) 



(a;,2/)P 



which involves a coupling realized via the entanglement transform K. Note that this forces an 
appropriate normalization of K, so as to guarantee the correct charge to flux ratio. 

3 Quantization of Hall resistance 



Note that if 



^1 

^2 



is a solution of (16), then in particular it is an eigenfunction of Ha- We now 



request that K he a stationary solution as in (|9|) with U = I, where 






are the eigenvectors of the Dirac Hamiltonian Ha- Assume u > E. Without loss of generality we 



may assume that 



^1 

i'2 



is on the list of mutually orthogonal eigenfunctions defining K via (9). 



(This can always be arranged, even if E has multiplicity greater than 1.) Thus, we have 

K 



i'2 



^2 



for a constant c. Consequently 



\K 




x,y)\'^ 



\H' + \H'- 



This implies that the last equation of ( 16 ) is equivalent to the statement 

d,A2-dyA,=RH{\H'+\i^2 



(17) 



Note that Rh is the ratio of the total flux to the total charge counted between the two types of 
particles. In particular, due to quantization of flux and charge, Rh is a rational number when 
expressed in appropriate units. We will demonstrate that in a symmetric realization of this model 
the Hall resistance is equal to 2Rh- We conduct the analysis in a rectangular-domain setting, and 
with the assumption that the wavefunctions are periodic in x, both with the same period, i.e. 



By substituting to (16), we readily obtain 



X2 



^1 


= exp 


{ikx)x 


y{y) 


^2 = 


= exp 


{ikx)x2 


(y)- 


uain 








= 


-ik 


-A,)x^ 


+ Ex2 


= 


{k- 


^l)X2 - 


-Exi. 



(18) 



Furthermore, let us assume A2 = 0, and Ai = Ai{y), Ai{0) = 0. Let us also introduce an auxiliary 
variable w = k — Ai. In such a case (16) is equivalent to the system 



w 

X'2 
By differentiating we readily obtain 



RHixl + X2), 

-wxi + Ex2 
WX2- Exi- 



w{Q) = k 



(19) 



x'l 

x'2' 



[-E^ + w^- w')xi 
[-E^ + w^ + w')x2- 



(20) 



These are 1-D stationary Schrodinger equations with energy E'^/2. Note that xi is affected by 
potential {w^ — w')/2, while X2 is affected by potential {w^ + w')/2. Both quantities result from 
the presence of the magnetic fleld, and thus represent the Hall effect. There are two potential 
functions, differing by the sign of the correction term w', because two types of carriers are involved 
in the charge transport. Note that 



w{y) = k + RHj{xl{y') + xliy'))dy'. 



(21) 



We calculate the current next. Recall that 

j = Reiijlii ■ dipi + ^lA)} + Re{4)l{i ■ #2 + ^2^)}- 
A direct calculation shows that the two l-forms amount to 

Re{ilj1{i ■ dilji + iPiA)} = -xKyMy)dx, 1 = 1,2. 

In summary 

1 1 w^ 

J = --^w'{y)w{y)dx = r— * d—. (22) 

(We remind the reader that * denotes the Hodge star.) Observe that the current flows along the 
X-axis. Since both types of carriers contribute to the current, the total Hall potential is the sum 
of the two separate potentials pointed out above, i.e. 

Vh = {w'^ + w')/2 + {w'^ -w')/2 = w'^. (23) 

It follows that 

J = ^ * dVu. (24) 

By integrating we then obtain that the Hall resistance is 2Rh- We remark that this is twice the 
value obtained in the Schrodinger type model, cf. |1], [5]. Recall that Rh has been introduced 
as the ratio of the number of quanta of the magnetic field to the number of quanta of electric 
charge e = 1 of both types of particles. If particles of type, say, ip2 were not accounted for in 
the calculation of the filling factor, then there would be a discrepancy between the inverse filling 
factor and the quantity Rh- 

It is interesting to ask when the electrons and holes are engaged as charge carriers in equal 
numbers. Let us position the edges of the plate symmetrically about the x-axis, i.e. y G [—h,h]. 
Note that the number of carriers of each type is proportional to 

xKy)dy\ 1 = 1,2. (25) 



Now assume X2(0) = Xi(0)) cind consider the last two equations of (19). In such a case one can 
find a solution satisfying the condition X2{,y) = Xi{~y)- ^oi this solution particles of each type 
occur in equal numbers. 

4 The elliptic system in the plane 

In the previous section we demonstrated that Hall resistance is quantized under certain symmetry 
assumptions, which reduce the problem to a system of ODEs. We emphasize that Hall potentials 



and Hall resistance can be defined and discussed in the general two-dimensional setting, [5]. In 
fact, the two-dimensional, i.e. PDE, setting cannot be avoided if we want to study the effect of 
lattice potential, etc. We will now briefly discuss the ramifications of the PDE problem. Consider 



(16) jointly with the assumptions about K put forth in the previous section. In this case it is 



equivalent to the following system of first order partial differential equations 

dyi/ji = id^il)i + {Ai + iA2)i)i + ^^2 

dytl,2 = -td,^2 - {Ai - tA2)^2 - Eiji (26) 

2 I L/, |2\ 



dyA, = d,A2-RH{\H'+\^: 



2| 



Note that the system will be strictly elliptic, if it is amended with a gauge condition 

dyA2 = -d,A,. (27) 



Note that the nonlinear terms in (26) are (real) analytic functions of the dependent variables. 
Therefore the classical Cauchy-Kovalevskaya Theorem, [U] , guaranties local existence of solutions 
of the initial value problem when the boundary conditions are also analytic. The initial value 
problem may help gain some insight into the nature of the system; an example of a solution 
obtained numerically with a finite difference method is displayed in Fig. 1. Since the boundary 
value problems spring to mind in this context, one should emphasize that in view of ellipticity, 
the natural setting is that of the Riemann-Hilbert type problem, [TD]. However, we wish to 
signal briefly that the boundary value problem is not the only type of approach appropriate for 
the study of the mesoscopic-loop model; there is a promising complementary approach. Namely, 



setting A = Ai+iA2, we can represent equations (26)-(27) as a system of three nonlinearly coupled 
9-equations: 

a^2 = -|(Av^2 + ^^i) (28) 

BA = -ii?^(|V;i|2 + |^2p). 
This suggests the relevance of a variety of complex methods, [TT], but this theme will not be 
discussed in the present article. 

Closing comments: The Dirac equation is undergoing a renaissance in condensed matter physics 
due to its relevance to the unusual two-dimensional crystal called graphene, e.g. [I2]. I emphasize 
that in this article we have not addressed any of the issues related to that stream of research. 
Rather, we have used the Dirac equation as an analytic alternative to the Schrodiner equation. It 
helped us find an almost equivalent reformulation of the nonrelativistic model. 
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Figure 1: A numerical solution of (26) with i? = 5 (in units adopted throughout the article), Rh 



5' 



and A2 = 0. The initial values '?/'2(a^,0) = 0.3exp(ix) + 0.2 * exp(2ia;), 'i/'i(a;, 0) = — 0.9-?/'2(a^,0), 



and Ai{x, 0) = are prescribed on the edge facing front-right. 



